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Abstract 

The meson spectroscopy is studied in the framework of the QCD string model. The mesons 
are composed of spinless quarks which move relativistically. We show that the QCD string model 
reproduces the straight line Regge trajectories and exibit the correct string slopes for both light- 
light and heavy-light mesons. We also present the exact numerical solution and the analytic 
approximation solution of the QCD string equations for the light-light mesons under the deep 
radial limit. We demonstrate that the QCD string model describes the universal relations for 
the light-light and heavy-light meson spectroscopies correctly from the comparison of the exact 
numerical solutions of the QCD string model for both light-light and heavy-light mesons. 



E-mail: henry ©physics . wise . edu 



1 



I. INTRODUCTION 



The mass spectrum of a hadron is a fundamental problem in hadronic physics. The 
QCD motivated potential model with scalar confinement in the nonrelativistic limit has 
been employed to deal with this problem. However the scalar potential confinement behaves 
badly 0] in the limit of small quark mass and diverges from QCD once we consider the 
relativistic corrections. It was found that the spin independent relativistic corrections using 
the low velocity Wilson loop formalism differed from those of scalar confinement |2j . The 
flux tube model which plays an important role in describing the relativistic meson states 
extends the nonrelativistic constituent quark model to include the gluonic degree of freedom. 
The review of the flux tube model may be found in Swanson's paper The QCD string or 
the relativistic flux tube(RFT) model which is the simplest generalization of the potential 
model overcomes these shortcomings of the scalar confinement. Its fundamental assumption 
is that the QCD dynamical ground state for large quark separation consists of a rigid straight 
tube-like color flux structure connecting the quarks and for massless quarks it reduces to 
Nambu string. This model is consistent with both spin-dependent and spin-independent 
QCD expectations while the scalar confinement potential model [5] is not consistent with 
spin-independent QCD relativistic correction. It includes the angular momentum and the 
rotational energy of the interacting field while only the field energy is considered in the 
potential model. Even though with this simple physical picture the RFT model coincides 
with the QCD relativistic correction derived from the rigorous Wilson loop formalism. Next 
we review the straight QCD string model and give the QCD string equations of mesons. 
Consider a light-light (LL) meson made up of two spinless quarks with equal mass mi. The 
quarks may rotate around their center of mass and vibrate in the radial direction. We obtain 
the following equations of the QCD string[6] for the LL mesons, 



2L 

— = 2W, n± v L + 2ar/K) (1) 
E = 2iy r7± + arS(v±) + V(r) (2) 

where I and E are angular momentum and energy of the meson respectively, V(r) rep- 
resents the short range potential, a is the confinement constant and r is the interquark 
separation. We also define the functions: 



arcsmv ± 

S(v ± ) = — — (3) 

fM = ^ r [S(v ± )-Jl~^] (4) 



W r = yjrf. + m 2 (5) 

Similarly, for a heavy-light meson(HL) made up of one heavy quark with mass M and 
one light quark with mass mi. They satisfy the conditions as M > m; -> and M>ar, 
we obtain the equations of QCD string for HL mesons, 
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- = W rl± v ± + 2arf{v ± ) (7) 
E = W r j ± + arS(v ± ) + V(r) (8) 

The paper is organized as follows. In section II we solve the RFT equations under two 
limits for the LL mesons. We demonstrate an analytical approximation solution under the 
deep radial limit in section III. The comparison of solutions for LL mesons with those for 
HL mesons is discussed in section IV. Special attentions are payed to the string slopes since 
they are significant criterions to tell whether a QCD motivated model is reasonable or not. 
Finally in section V our conclusions are summarized. 



II. SOLUTIONS UNDER TWO LIMITS 
A. Circular Motion limit 

We now consider the circular motion limit for the LL mesons in which the two light 
spinless quarks moving circularly around their center of mass. So this will reduce the p r — > 
and v± — *■ 1. This reduction yields 

SM = I (9) 
/K) = \ (io) 

Pluging functions S(v±), f{v±) back into the general LL RFT equations, we obtain 

2_L = vror 

r 4 
nar 

E = 12 

2 V ; 

Therefore we find the 'Narnbu' slope for the LL meson is, 

Q " = = 2^ (13) 
We note that it is equal to half of that for HL meson. We will discuss this slope again 
in section IIVI where the results can be seen clearly in the corresponding graphs no matter 
whether there is color Coulomb potential or not. 



B. The deep radial limit solutions for LL mesons 

Under the deep radial limit, considering the LL mesons and assuming again the light 
quark mass is zero we want to explore the effect of the short range color Coulomb potential. 
That the QCD string reduces to a time component vector potential was shown under the 
deep radial limit in which L <C i? 2 .|7] It is reasonable to apply the semi-classical quantization 
condition for a spherically symmetric system to quantize the QCD string, 
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1 f r + 

(n + -)tt = / drp r = I (14) 

Z J V — 

where r + ,r_ represent the radial distances at the turning points of the motion. Assume 
the potential takes the simple form V(r) = —-[8]. we obtain under the deep radial limit, 



4p 2 r + A^ = (E-ar + ^, (15) 



p r = —^{Er - ar 2 -2L + k)(Er - ar 2 + 2L + k). (16) 



We now define the dimensionless variables, 



x = \far (17) 

(18) 

(19) 



2^/E 

Pr 

fa 



Then we obtain using the new dimensionless variables 
P 



' " ' 'a 2x 2.r '' 



where 



Q 1 = x 2 - 2ex + 2L - k, (21) 
Q 2 = x 2 - 2ex — 2L — k. (22) 



The turning points follows: 



x x = e - Ve 2 + 2L + k, (23) 
x 2 = e - Ve 2 -2L + k, (24) 
x 3 = e + Ve 2 -2L + k, (25) 
x 4 = e + Ve 2 + 2L + k, (26) 

(27) 



Now we are ready to evaluate the integral 



/ S r drp, .= r dx^^ (28) 

Jr. Jx 2 2x 

Assume, 

x = uy + v (29) 

therefore, 
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Qi = {uy+v—X2){uy+v — xz) = u 2 y 2 +uy{2v — X2 — x 3 ) + {v—X2){v — xz) = u 2 y 2 — (e 2 — 2L+k) 

(30) 



where we let 



x 2 + x 3 



(31) 



Similarly, 



Q2 = (uy+v — Xi)(uy+v— X4) = u y +uy(2v—x 1 —X4) + (v—xi)(v—X4) = u y — (e +2L+k) 

(32) 



Define 



then, 



u 2 = e 2 - 2L + k 



(33) 



Q1Q2 = u 2 (l - y 2 )(l - my 2 )(e 2 + 2L + k) 



(34) 



where 



Due to the following relation, 



m 



e 2 -2L + k 
e 2 + 2L + k 



(35) 



x 2 — e = —u 
x 3 — e = u 



(36) 
(37) 



plug them into Eq. 



Jx 2 2x 
i 2 r 1 dy 



u 



2 J-i uy + e 

iiVe 2 + 2L + k r fO dy 



(1 -y 2 ){l -my){e 2 + 2L + k) 



1 uy + e 



(1 — y 2 )(l — my 2 ) + 



dy 



uy + e 



(1 - y 2 )(l - my 2 ) 



u 2 Ve 2 + 2L + k r 1 



d y y(l-y 2 )(l-my 2 )( ^— + —L- ] 

L v — uy + e uy + e 



(38) 



Therefore we get the TCV integral for the LL mesons under the deep radial limit, 
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_ (e 2 -2L + k)Ve 2 + 2L + k r 1 
e 7o 



'1 7/2V1 e - 2L + fc 2 > 
K e 2 + 2L + fc y J 

l- e2 -2f + V 



( e 2_2L + fcye 2 + 2L + fc j 



where 



J = / 2 9r 7 r • (40) 

Jo i _ e — 2L + k 2 

e 2 y 

Using the semi-classical quantization condition for a spherically symmetric system to 
quantize the QCD string, we obtain, 



'(i-z/)(i- %T^ + ^ 



(e 2 -2L + fc) V / (e 2 + 2L + A;) ,i V ll " ^ ni "?T2LTI y ^ , 1, (| 

/ = / dy~ 2 — or , 7 = {n + -)^ 41 

e Jo 1 e - 2L + fc „, 2 2 

1 ? y 

When evaluating Eq. (j41|) we need replace the classical angular momentum with Langer 
correction to take into account the centrifugal singularity, namely L — > Moreover, this 

integral I may be represented with complete standard elliptic integrals, 



(e 2 — 2L + k)J (e 2 + 2L + k) ml ml 1 

I = K - ^ L[(l - -)(1 - -)U(h, m) + -(1- -)F(m) + -E{m)} (42) 

where 

e 2 -2L + A; 

/i= ^ . (43) 

III. ANALYTIC APPROXIMATION SOLUTION 

Next we derive an analytical approximation solution of QCD string equations for LL 
mesons under the deep radial limit L <C E 2 . After some tedious calculations (Shown in Ap- 
pendix) we obtain an analytical approximation of the spectroscopic relation for LL mesons, 



We now compare this approximation result with the exact numerical solution Eq. (j41jl . 
We draw the 6 figures for every other n starting from to 10. From these figures(from 
FIG.l to FIG. 6) we easily identify that the approximations are more accurate along with 
the increasement of n since the deep radial limit is better satisfied. 
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FIG. 1: ground state n=0 



IV. COMPARISON OF SOLUTIONS FOR LL MESONS WITH THOSE FOR HL 
MESONS UNDER THE DEEP RADIAL LIMIT 

Next we take advantage of numerical methods to draw the graphs for LL meson integral. 
At first we turn off the color Coulomb potential. For the ground state, the 1st and 2nd 
excited states, the relation of angular momentum and energy is shown respectively in FIG|7| 

From the FIGEIwe easily notice that the slopes of the curve for LL meson are half of 
those for HL meson which are shown in FIG|H1 The RFT mode reproduces the straight line 
Regge trajectories and exibit the correct string slopes for both LL and HL mesons. 

Now considering the effect of color Coulomb potential and assuming that k = 0, k = 0.2 
and k — 0.4 separately, then we obtain the curves modified by k values in FIG El Again we 
find that the slopes of the curve for LL mesons are half of those for HL mesons shown in 
FIGHU1 

V. SUMMARY AND CONCLUSION 

We have studied the relativistic flux tube model and its applications to the meson spec- 
troscopies. Under two limit we also solve the relativistic flux tube equations for LL mesons. 
We have presented the exact numerical solution for light-light mesons and compared the 
numerical solution with an analytical approximation solution in the deep radial limit. The 
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FIG. 2: The 2nd excited state 



comparison of solutions for LL mesons with those for HL mesons is also discussed. 

We conclude that the RFT model reproduces the straight line Regge trajectories and 
exibit the correct string slopes for both LL and HL mesons. They appear very naturally in 
the RFT approach. The RFT model correctly account for the rotation of the QCD string. 
From the comparison of exact numerical solution with the analytic approximation solution 
for light-light mesons, we find that they are much more in agreement with each other for 
higher excited states since the deep radial limit is better satisfied. The relativistic flux tube 
model describes the universal relations for LL and HL meson spectroscopies correctly. 
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Appendix Derivation of the analytical approximation solution under deep radial 
limit 

We now derive an analytic approximation solution of RFT model for LL mesons under 
the deep radial limit L <C E 2 and k — 0. To this purpose, we define 
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FIG. 3: The 4th excited state 
Since [3 <C 1 in the deep radial limit, m and h may be represented with f3 as 



in --= - — £ « 1 — 2/3 



1 + /3 



Plugging into Eq.(|4*0]) we obtain 



-1 y/(l-y*)(l-y2 + 2l3y*) 

J = I ay- ~2 — —2 



o 



Assuming that, 



i-y 2 



l 



1 + X 



then J may be represented with the new parameters x and f3 as 



1 f°° 
J = - dx- 

2 Jo 



1 + 2(3x) 



\ x(l + 



x 



(46) 
(47) 



(48) 
(49) 



l + (50) 
Next expand J with taylor expansion and drop terms of second order in f3 or higher, so 
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FIG. 4: The 6th excited state 



l r°° l 
J(0) = 2i 



(i + xY 



(51) 



(52) 



dJ_ 
~dp 



(3=0 



1 r dx 1 

2 Jo V x(l + x) 3 



1 + /3a; ; 



dp 



dz- 



1 



1 7T 



2^o y(TT2^(l + ^) 2 2 4 
where we let z = fix. Therefore when < 1, we find to the first order in P 



(53) 



_ (e 2 -2L)Ve 2 + 2L 



J = e *(l-/?)(l + /J)*[l + (± 



n 



= ^(1 - t/5) 



(54) 



After substituting e 2 with we obtain an analytical approximation of the heavy-light 
semi-classical time component vector potential integral, 
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FIG. 5: The 8th excited state 





^<l-;« 

2a 4 J 



(55) 



Applying the semi-classical quantization condition for a spherically symmetric system 
and replacing the classical angular momentum with Langer correction to take into account 
the centrifugal singularity, we immediately obtain the spectroscopic relation 



L + 2n + - 
2 



E 2 
2na 



(56) 



which is the Eq. ()44j) and an excellent approximation to the numerical solution of RFT 
model for a LL meson presented in section IIHI 
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FIG. 6: The 10th excited state 
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FIG. 7: The middle dashed line is for the first excited state while the right one is for the 2nd 
excited state. They compare with the ground state graph represented by the solid line. 
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FIG. 8: The middle dashed line is for the first excited state, the right one is for the 2nd excited 
state while the solid line is for the ground state. 
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FIG. 9: The left dashed curve corresponds to the color Coulomb constant k = 0.4, the middle 
dashed curve is for k = 0.2 while the solid one is for k = 0. 
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